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In some processes at the LHC, theoretical precisions of 1% are desired. With an eye toward such 
precisions, we introduce the theory of the simultaneous YFS resummation of QED and QCD to 
compute the size of the expected resummed soft radiative threshold effects in precision studies of 
heavy particle production at the LHC. Our results, that the soft QED threshold effects are at the 
level of 0.3% whereas the soft QCD threshold effects enter at the level of 20%, show that both must 
be controlled to be on the conservative side to achieve such precision goals. 



1 Introduction 

In high energy collider processes, such as tt 
production at FNAL, precision predictions 
for soft multiple gluon (n(g)) effects are al- 
ready needed: the uncertainty on m t [1], 
8mt = 4.3 GeV, receives a soft n(g) un- 
certainty ~ 2-3 GeV, for example. At the 
LHC/ILC, the requirements will be even 
more demanding and soft n(g) MC exponen- 
tiation results will be an important part of 
the necessary theory - YFS exponentiated 
0{a 2 s )L calculations, in the presence of par- 
ton showers, on an event-by-event basis. 

How relevant are QED higher order cor- 
rections when QCD is controlled at ~ 1% 
precision? Many authors [2] are preparing 
the necessary results that would lead to such 
a precision on QCD for LHC processes. Es- 
timates by Refs. [3-7] show that one gets 
few per mille effects from QED corrections 
to structure function evolution. The well- 
known possible enhancement of QED cor- 
rections at threshold, especially in resonance 



production, leads us to estimate how big are 
these effects at the LHC. 

We treat QED and QCD simultaneously 
in the respective YFS [8,9] exponentiation to 
estimate the role of the QED threshold ef- 
fects at the LHC in the representative pro- 
cesses pp — > V + 71(7) + m{g) + X — ► It' + 
n'(-f) + m(g) + X, where V = W ± ,Z,and 
I = e,n, £' = v e , Vp{e, /x) respectively for 
V = W+{Z), and £ = v e ,v^, £' = e,/i re- 
spectively for V = W~ . Precision studies of 
these processes have been proposed for lumi- 
nomctry at the LHC [10] and at FNAL [11], 
where 2-3% is the target precision tag for the 
LHC, for example. The latter would indeed 
require a theoretical precision tag of ~ 1% 
in order that the theory error not figure too 
prominently in the over-all precision. 

Our discussion is organized as follows. 
After giving a brief review of the YFS theory 
and its extension to QCD in the next sec- 
tion, in Section 3 we introduce QED®QCD 
YFS exponentiation. In Section 4, we apply 
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the new development to the threshold correc- 
tions in single V production at the LHC and 
at FNAL. Section 5 contains some concluding 
remarks. 

2 Review of the YFS Theory and 
its Extension to QCD 

As realized in Refs. [9] by Monte Carlo 
methods, for e+(pi)e~(<?i) ->• /(p2)/(ga) + 
71(7) (fci, •, k n ), renormalization group im- 
proved YFS theory [12] gives, 



da 



exp 



„2a ReB+laB v^°o 



n=0 ■' 11.7=1 ~~kT 

f d 4 y iy(pi+qi-p2-q2~T. a kj)+D 



, ^ a p 2 d q-2 /i\ 



where the YFS infrared functions B, B and 
D are known. For example, the YFS hard 
photon residuals (3i in (1), i = 0,1,2, are 
given in the first paper in Ref. [9] and re- 
alize the YFS exponentiated exact 0(a) and 
LL 0(a 2 ) cross section for Bhabha scattering 
via a corresponding Monte Carlo realization 
of (1). 

In Refs. [13, 14] we have extended the 
YFS theory to QCD: 



da c 



n _ SUM m (QCD) 



n=0 j=l 



d 3 k, 



! 



d V iy(Pi+P2-Qi-Q2-T, fcj)+£>QCD 
(2^)4 



d i P 2 d A Q 2 
Q2 



(2) 



where gluon residuals n (ki, . ■ ■ , k n ) , defined 
by Ref. [13], are free of all infrared diver- 
gences to all orders in a s (Q). The functions 
SUM IR (QCD), D QCD , together with the ba- 
sic infrared functions Bq^ d , Bqc D , Sqcd 
are specified in Ref. [13]. We call atten- 
tion to the essential compensation between 
the left over genuine non-Abelian IR virtual 
and real singularities between J dPh(3 n and 
J dPh(3 n +i respectively that really allows us 



to isolate (3j and distinguishes QCD from 
QED, where no such compensation occurs. 

We stress that the YFS resummation 
which we exhibit here is fully consistent with 
that of Refs. [15, 16]. We refer the reader to 
Ref. [17] for more discussion of this point. 

3 Extension to QED<g>QCD and 
QCED 

Simultaneous exponentiation of QED and 
QCD higher order effects [17] gives 



QCD 



B 
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QCED, 
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which leads to 



d<Toxp = e 



SUM IR (QCED) 
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(4) 



where the new_ YFS residuals, defined 
in Ref. [17], /3 nm (ki,...,k n ;k' 1 ,...,k' m ), 
with n hard gluons and m hard pho- 
tons, represent the successive application 
of the YFS expansion first for QCD 
and subsequently for QED. The functions 
SUMir(QCED), -Dqced are determined from 
their analoga SUMir(QCD), £>qcd via the 
substitutions in (3) everywhere in expressions 
for the latter functions given in Refs. [13]. 

Infrared Algebra(QCED): the average 
Bjorkcn x values 

x aV g{QED) = j(QED)/(l + j(QED)) 
x avg (QCD) = 7 (QOD)/(l + ~f(QCD)) 

where 7(A) = *±£±(L a - 1), A = 
QED, QCD, with Ca = Q% Cf, respectively, 
for A = QED, QCD and the big log L s , im- 
ply that QCD dominant corrections happen 
an order of magnitude earlier than those for 
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QED.This means that that the leading /3q°q°' ) - 
level gives a good estimate of the size of the 
effects we study. 

4 QED®QCD Threshold 
Corrections at the LHC 

We shall apply the new simultaneous 
QED®QCD exponentiation calculus to the 
single Z production with leptonic decay at 
the LHC ( and at FNAL) to focus on the 
1SR alone, for dcfiniteness. See also the work 
of Refs. [18-20] for exact 0(a) results and 
Refs. [21-23] for exact 0(a 2 s ) results. 

For the basic formula (we use the stan- 
dard notation here [17]) 



da eX p(pp -> V + X ->&' + X 1 ) 



(5) 



we use the result in (4) here with semi- 
analytical methods and structure functions 
from Ref. [24]. A Monte Carlo realization 
will appear elsewhere [25]. 

We do not attempt to replace HER- 
WIG [26] and/or PYTHIA [27] - we intend to 
combine our exact YFS calculus with HER- 
WIG and/or PYTHIA by using the latter 
in lieu of the {Fi}. This combination of 
theoretical constructs can be systematically 
improved with exact results order-by-ordcr 
in a s , where currently the state of the art 
in such a calculation is the work of Frixionc 
and Webber in Ref. [28] which accomplishes 
the combination of an exact 0(a s ) correction 
with HERWIG. We note that, even in this 
latter result, the gluon azimuthal angle is av- 
eraged in the combination. We note that the 
recent alternative parton shower algorithm 
by Jadach and Skrzypek in Ref. [29] can also 
be used in our theoretical construction here. 
Due to its lack of the appropriate color co- 
herence [30], we do not consider IS A JET [31] 
here. 

We compute , with and without QED, 
the ratio : 



i exp 




(We stress that we do not use the narrow res- 
onance approximation here.) 



, QCED = QCD+QED, LHC 
,QCD, LHC 

, QCED = QCD+QED, Tevatron 
,QCD, Tevatron. 

(6) 

We see that QED is at the level of .3% at 
both LHC and FNAL. This is stable under 
scale variations [17]. We agree with the re- 
sults in Refs. [18-22] on both of the respective 
sizes of the QED and QCD effects. The QED 
effect is similar in size to structure function 
results in Refs. [3-7]. 



exp 



o~expl o~Bom to get the results 



5 Conclusions 

YFS theory (EEX and CEEX) extends to 
non-Abelian gauge theory and allows simul- 
taneous exponentiation of QED and QCD. 
For QED(g)QCD we find that full MC event 
generator realization is possible in a way 
that combines our calculus with Hcrwig and 
Pythia in principle. Semi-analytical results 
for QED (and QCD) threshold effects agree 
with literature on Z production. As QED is 
at the .3% level, it is needed for 1% LHC 
theory predictions. A firm basis for the com- 
plete 0(a 2 , aa s , a 2 ) results needed for the 
FNAL/LHC/RHIC/TESLA/ILC physics has 
been demonstrated and all of the latter are in 
progress. 
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